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Abstract. We discuss the potential scattering on the noncompact star graph. The 
£S| ' Schrodinger operator with the short-range potential localizing in a neighborhood of 

the graph vertex is considered. We study the asymptotic behavior the corresponding 

H . ' scattering matrix in the zero-range limit. It has been known for a long time that in 

{/*} • dimension 1 there is no non-trivial Hamiltonian with the distributional potential 5', 

fH | i.e., the 5' potential acts as a totally reflecting wall. Several authors have, in recent 

years, studied the scattering properties of the regularizing potentials ae~ 2 Q(x/e) 
approximating the first derivative of the Dirac delta function. A non-zero transmission 
through the regularized potential has been shown to exist as e — > 0. We extend these 
results to star graphs with the point interaction, which is an analogue of 5' potential 
on the line. We prove that generically such a potential on the graph is opaque. We 
qq ' also show that there exists a countable set of resonant intensities for which a partial 

0^ ■ transmission through the potential occurs. This set of resonances is referred to as 

ilj | the resonant set and is determined as the spectrum of an auxiliary Sturm-Liouville 

problem associated with Q on the graph. 
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1. Introduction 



Ctf 



Schrodinger operators with potentials supported on a discrete set of points (such 
potentials are usually termed "point interactions" ) have attracted considerable attention 
both in the physical and mathematical literature over several past decades. One of 
the reasons for this is that such singular Hamiltonians are widely used in various 
application to atomic, nuclear, and solid state physics. Applications also arise in optics, 
for instance, in dielectric media where electromagnetic waves scatter at boundaries or 
thin layers. Another reason is that Schrodinger operators with point interactions often 
form "solvable" models in the sense that the resolvents of such operators can explicitly 
be calculated. Consequently the spectrum, the eigenf unctions, as well as resonances and 
scattering quantities, can also be determined explicitly (see [2j[3] and references therein). 
In the physically oriented literature point interactions are often understood as sharply 
localized potentials, exhibiting a number of interesting features and unusual effects not 
seen for regular potentials. 



On 8' -like potential scattering 2 

Currently, there is increasing interest in solvable models on graphs in particular, as a 
reaction to a great deal of progress in fabricating graph-like structures of a semiconductor 
material, for which graph Hamiltonians represent a natural model (see the survey [20] 
for details). The idea to investigate quantum mechanics of particles confined to a graph 
originated with the study of free electron models of organic molecules [2H - I26] . Among 
the systems that were successfully modeled by quantum graphs we mention e.g., single- 
mode acoustic and electro- magnetic waveguide networks [H], Anderson transition [I] 
and quantum Hall systems with long range potential (6], fracton excitations in fractal 
structures [5], and mesoscopic quantum systems pj3]. The essential component of such 
graph models is the wavefunction coupling in the vertices. The interface conditions 
have to be chosen to make the Hamiltonian self-adjoint, or in physical terms, to ensure 
conservation of the probability current at the vertex. 

One of the most natural way to define a graph Hamiltonian corresponding to a point 
interaction supported by the branching point is to choose an appropriate self-adjoint 
extension of the corresponding free Hamiltonian with the interactions points removed. 
This approach has been realized in [13], where the authors found that a vertex joining 
n graph edges can be described by n 2 real parameters defining the interface condition 
at the vertex. A general form of such a coupling was described in [18] by a pair of nxn 
matrices A, B such that rank (A, B) = n and AB* is self-adjoint; the boundary values 
have to satisfy the conditions 

A$(a) + B*'(a) = 0. (1) 

Here the symbol ^(a) is used for the column vector of the boundary values at the vertex 
a, and analogously \l/'(a) stands for the vector of the first derivatives, taken all in the 
outgoing direction. 

Another possible approach to define the vertex coupling is to approximate a 
quantum graph with a thin quantum waveguide (see [H[8l[T2l|23] and references therein). 
Different vertex couplings were also discovered in [71 fT0llTT| ri3]. 

The Schrodinger operators with the Dirac delta-function and its derivatives in 
potentials have been studying intensively since the eighties of last century (see [2~|l3jl9] and 
the references given there). One of the first well-studied Hamiltonian was the operator 
with Dirac's delta-function potential 5. A case of special interest arises when the 5- 
potential is replaced by the derivative 5' of the Dirac delta-function. Seba [27] appears 
to have been among the first to consider such a Hamiltonian. To define the Hamiltonian 
— ^2 + 5'(x) he approximated 5' by regular potentials e~ 2 Q{x/e) and then investigated 
the convergence of the corresponding family of regular Schrodinger operators 

d 2 1 „/x\ 



dx 2 e 2 \ e . 

Here Q is an integrable function. If Q has zero mean and the integral j R t Q(t) dt is equal 
to —1, then the sequence e~ 2 Q(x/e) converges in the sense of distributions as e — > to 5', 
which served as motivation for such an approximation. Seba claimed that the operators 
S £ converge in the uniform resolvent sense to the direct sum of the unperturbed half-line 
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Schrodinger operators subject to the Dirichlet boundary conditions at x = 0. From the 
viewpoint of scattering theory it means that the <5'-barrier is completely opaque, i.e., 
in the limit e — > the potential e~ 2 Q(x/e) becomes a totally reflecting wall at x = 
splitting the system into two independent subsystems lying on the half-lines (— oo, 0) 
and (0, oo). 

But this result is in conflict with conclusions reached in [28], where the resonances 
in the transmission probability for 5'-like potential have been observed. In that paper 
an exactly solvable model with a step function aQ was considered. The authors found 
a discrete set of intensities a n , for which partial transmission through the limiting 8'- 
potential occurs. The values a n are roots of a transcendent equation depending on the 
regularization Q. Exactly solvable models with other piecewise constant potentials as 
well as nonrectangular regular izat ions of 8' have later been studied in [29H32] and the 
same conclusion has been drawn. 

In [T5l[T6] the authors approximated the formal Hamiltonian — J^ + q(x) + a8'(x) 
on the line by regular Schrodinger operators 

d 2 



S £ (a, Q) = -— + q(x) + ^q(^ 



and the similar resonant effect was discovered. Here the function Q G C^°( — 1, 1) has 
zero mean, a is a real coupling constant, and q is a real valued fixed potential tending 
to +oo as |x| — > oo, which ensures that the spectrum of S £ (a,Q) is discrete. The map 
assigning a selfadjoint operator S(a,Q) to each pair (a, Q) was constructed there. The 
choice of S(a,Q) was determined by proximity of its energy levels and pure states to 
those for the Hamiltonian with regularized potentials for small e. It was established 
that for almost all coupling constants the operator S(a, Q) is just the direct sum of the 
Schrodinger operators with potential q on half-axes subject to the Dirichlet boundary 
conditions at the origin. But in the exceptional case the nontrivial coupling at the origin 
arises. The notion of the resonant set Eq, which is the spectrum of the problem 

-g" + aQ(t)g = 0, te (-1, 1), g'(-l) = </(l) = 0, 

with respect to the spectral parameter a was introduced in [13 (TB]. It was also 
constructed the coupling function 9q : Sq — > M. defined via #q(«) = g a (l)/g a (— 1), 
where g a is an eigenfunction corresponding to the eigenvalue a G Sq. In the case 
when the coupling constant a belongs to the resonance set, the ^'-barrier admits a 
partial transmission of a particle, and an appropriate wave function obeys the interface 
conditions 

V>(+o) = m«M-o), e Q (aW(+o) = ^(-o). 

Studies of [151 [16] have been continued in [T7J[2T]. First the findings of [28H32] 
were generalized in [21], where the scattering on an arbitrary potential of the form 
ae~ 2 Q(x/e) has been considered. It was proved that such a potential is asymptotically 
transparent only if a coupling constant belongs to the resonant set. Moreover, the 
scattering amplitude for the Hamiltonians S s (a,Q) and — ^ converges as e — > to 
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that for the limiting Hamiltonian S(a,Q). In [17] the authors have not only pointed 
out a mistake in [27], but also have showed that the operators S e (a, Q) converge in the 
uniform resolvent sense as e — > to S(a,Q). The results of [T71I2T] were obtained for 
the special case q = 0; generally, the same can be derived without difficulty 

Many of point interactions on the line were extended to deal with graphs. One of 
the best known is the 8 coupling at n edge vertex: 

n 

V>i(a) = . . . = ip n (a), ^2 ^i(°) = a ^(°) 

(see, e.g., pT0 | [TT]). Such a model is a generalization of the Hamiltonian 

given by S a> sf = —f" on the domain 

V(S a , s ) = {fe H 2 (R \ {0}) I /(+0) = /(-0), /'(+0) - f'(-0) = o/(0)}. 

A similar generalization is possible for the model for Scrodinger operators with 8'- 
interactions, defined by Sp#f = —f" on the set of functinctions 

V{S Pfi ) = {feH 2 (R\ {0}) | /' (+0) = f(-0), /(+0) - /(-0) = fif (0)}. 

This generalization was described in [T0|fTT]: 

n 

^(a) = ... = <(a), ^^(«) = W(«), 

i=i 

and was called 8 ; s coupling. In the following section we briefly sketch the findings of [22], 
where an analogue for the Hamiltonian S e (a,Q) was considered on the metric graph. 

1.1. Graph Hamiltonian with the 8' -like potential 

It will be convenient to recall basic notions of the theory of differential equations 
on graphs. By a metric graph G = (V, E) we mean a finite set V of points in R 3 
(whose elements are called vertices) together with a set E of smooth regular curves 
connecting the vertices (the elements of E are called edges). The sets of vertices and 
edges of the graph G are sometimes denoted by V(G) and E(G) respectively. A map 
/ : G —$■ R is said to be a function on the graph, and the restriction of / on the edge 
g G E{G) will be denoted by f g . Let us introduce the space of functions on the graph 
C°°(G) — {/ | f g e C°°(^) for all g G E(G)}. Each edge is equipped with a natural 
parametrization; the differentiation is performing with respect to the natural parameter. 
We denote by df/dg(a) the limit value of the derivative at the point a G V(G) taken 
in the direction away from the vertex. The integral of / over G is the sum of integrals 
over all edges 



I fdG= J2 [f d 9- 



g€E(G) 
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Let L 2 (G) be a Hilbert space with the inner product (u,v) = J G uvdG and with 
the norm || • ||l 2 (g)- We also introduce the Sobolev spaces H k {G) = {/ G L 2 (G) \ 
f g G H k (g) for all g G E(G)} and the space of absolutely continuous functions 
AC(G) = {f\f g e AC(g) for all g G E(G)}. 

Let us consider a noncompact star graph T consisting of three edges 71, 72 and 
73. All edges are connected at the vertex a. Then E(T) = {71,72,73} and suppose 
that all edges are half-lines. We write ah for the point of intersection of the e-sphere, 
centered at a, with the edge jj G E(T) and denote by T e a sub-partition of T containing 
new vertices a\, a\ and a\. Each a 6 , divides the edge jj of the graph T into two edges 
Uj and 7J of the graph T e (see figure [p. Let Q £ be a star subgraph of T e such that 
V(Q £ ) = {a,af,a|,a|} and E(Q £ ) = {o;f,a;|,a;|}. 




Ti 



Yi 




Q 



CO, 



Figure 1. Graphs T e and Q 



By fi C It we denote the e _1 -homothety of the graph Q £ , centered at a. Obviously, 
the resulting graph does not depend on the small parameter e. The graph Q is a star 
with a center at the origin b of the auxiliary space M? and with the vertexes a,j, which 
are the images of the points ah. The edge of fl connecting b to a,j will be denoted by uij. 

We introduce the set Q = {Q G C°°(Q) \ f n Qdfi = 0}. For each nonzero element 
Q G Q let us define the sequence 

:~ 2 Q((x — a)/e) for x G Q £ , 

for x G T £ \ Q £ . 

The potential Q £ is referred to as the 5'-like potential. 

In [22] , it was discovered the family of Schrodinger operators on the star graph 

R £ (a,Q)f = -f" + (q + aQ £ )f, 

V(R £ (a, Q)) = {fe L 2 (T) \ f, f G AC(T), -f" + (q + aQ £ )f G L 2 (r), 



Q £ {x) 



Li(a) = /o,|(a) = / w «(o), 



E^< 



1 duj ) 







where Q e Q, a is a real coupling constant, and q is a smooth real- valued potential 
such that q y (x) — > +00 as |x| — > +00 for all 7 G -E(r). Such behavior of g ensures 
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the discreteness of the spectrum of H e (a, Q). A self-adjoint operator H(a, Q) has been 
assigned to each pair (a, Q). The choice of an operator is argued by the close proximity 
of the energy levels and the pure states for the Hamiltonians H e (a,Q) and H.(a,Q) 
respectively. Two spectral characteristics of the function Q are introduced in [22] : the 
resonant set S, which is the spectrum of the eigenvalue problem 

-g" + aQg = on Q\V(Q), (2) 

3 , 
&*(&) = 9oj 2 (b) = gu 3 (b), Yl d^~( b } = °' ( 3 ) 

dg , dg dg 

d^ {ai) = A^ (a2) = A^ (a3) = °' (4) 

and the coupling function : S — > CP 2 , where CP 2 is the complex projective plane. The 
spectrum of the problem (T5]) is real, discrete and has two accumulation points ±oo. It 
consists of simple and double eigenvalues. Write Si for the subset of simple eigenvalues 
of the problem ([2]) and by S 2 we denote the subset of double eigenvalues. The resonant 
set can be represented as the union of Si and S 2 . 

Let a be a simple eigenvalue of the problem (T5]) with an eigenfunction u a such that 
u a( a i) + n a( a 2) + u l( a 3) = 1, then introduce 

0i (a) = iia(ai), 2 (a) =u a (a 2 ), 9 3 (a) = u a (a 3 ) . 

For a G S 2 we consider 

0i (a) = f Q (a2)^a(a3) - v a (a 3 )w a (a 2 ), 
02(a) = ^a(>3)w a 0i) - w a (ai)w a (a 3 ), 
03(a) = ^a(ai)w a (a 2 ) - v a (a 2 )w a (a 1 ), 

where v a and w a form a base in the corresponding eigenspace such that 2 (a) + 0|(a) + 
0|(«) = 1. The coupling function is defined via 0(a) = (0i(a), 2 (a), 03(a)). It if easy 
to check that if we change a base in the eigenspace the point (0i(a), 2 (a), 03 (a)) G C 3 
should be replaced by (A0i(a), A0 2 (a), A03(a)) with some complex A, hence we shall 
find it convenient to consider as a function from the resonant set into CP 2 . In the 
non-resonant case, when the coupling constant a does not belong to the resonant set, 
H(a, Q) is the direct sum of the Schrodinger operators with the potential q on edges, 
subject to the Dirichlet boundary conditions at the vertex a. If a G Si (the case of 
simple resonance), then the operator H(a,Q) acts via H(a, Q)ifj = —ip" + qtp on an 
appropriate set of functions obeying the interface conditions 

2 (a)0 3 (a)^ 71 (o) = 0i(a)0 3 (a)^ 72 (a) = 0i(a)0 2 (a)^ 73 (a), 

0i (a) %{a) + 2 (a) J|(a) + 3 (a) %(a) = ^ 

that can be expressed by (CQ) with 

/ 3 (a) -0i(a) \ / \ 

A= 3 (a) -82(a) , B 

^0 j 





y 0i(a) 2 (a) 3 (a) J 
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In the case of double resonance (a G S 2 ), the interface conditions may be written as 

6 2 (a)6 3 (a)^(a) = e l (a)e 3 (a)^(a) = 0i(a)fc(a)|£(o), 
9i(a)ip 71 (a) + 9 2 (a)xjj l2 (a) + 3 (a)ip l3 (a) = 

To derive these conditions in the form ([1]) we need only interchange the roles of the 
above matrixes A and B. 

If A is an eigenvalue of the operator H(a, Q), we will denote by Pa the orthogonal 
projector onto the corresponding eigenspace. Let P\(s) stand for the orthogonal 
projector onto the finite dimensional space spanned by all eigenf unctions corresponding 
to those eigenvalues \ £ of H £ (a, Q) that A e — > A as e — > 0. The results of [22] may be 
summed up in the following theorem. 

Theorem 1.1 All eigenvalues ofH £ (a,Q) (except at most a finite number) are bounded 
as e — > 0. Let X £ be an eigenvalue of H. £ (a,Q) bounded as e -> 0, then X £ has a finite 
limit A that is a point of the spectrum ofH(a,Q). Moreover, 1 1 P_x (e) — Pa || — >• as 
e — > 0. Conversely, if X is an eigenvalue ofH(a,Q), then there exists an eigenvalue X £ 
of H e (a, Q) such that X £ — > X as e — > 0. 

Although it has been showed the close proximity of the energy levels and pure 
states for the limiting and regularized Hamiltonians, we are still in the dark about 
convergence as e — > of the operators H. £ (a,Q) in any topology As a result, we do 
not know anything about convergence of the scattering quantities or other physical 
characteristics. Our objective in this paper is to give one further motivation for the 
choice of the Hamiltonian H(a,Q). We shall study the scattering properties of the 
finite-range potentials aQ £ on the graph T £ in the limit e — > 0. We prove that the 
scattering coefficients depend on the intensity a and the function Q in such a way that 
for all values of a the barrier aQ £ is completely opaque except for the set S of resonant 
values, at which a partial transmission through the potential occurs. It will also be 
shown that the scattering amplitude for the Hamiltonians H £ (a,Q) and Ho converges 
as e — > to that for the limiting Hamiltonian H(a,Q). Here the Hamiltonian Ho of a 
free particle on T acts via H / = — /" on its domain consisting of those functions from 
H 2 (T) that are continuous on T and satisfy the Kirchhoff boundary conditions at the 
vertex. 

2. Scattering problem for the <5'-like potential on the graph 

Let us state the main result of this paper. 

Theorem 2.1 For each k > and a 6 1 the scattering matrix for the operators 
H e (a, Q) and Ho converge as e — > to the scattering matrix for H(«, Q) and Ho. 

To start with, we briefly treat stationary scattering on the graph T associated with 
the Hamiltonians H(a, Q) and Ho- From now on we assume that q is a zero function, i.e., 
the operator H(a;, Q) involves no potential and H £ (a, Q) has only perturbed potential 
aQ £ . It is sufficient for us to look at the nontrivial case when the coupling constant a 
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belongs to the resonant set £ = Si U S 2 , since in the opposite case the particle is always 
reflected. Let us introduce a natural parametrization s G (0, +oo) on each edge of T and 
T e , where s = corresponds to the vertex a. Consider the incoming monochromatic 
wave e~ lks with k > coming from infinity along the edge j n . The corresponding 
wavefunction has the form 

T nm e lks on j m , m^n, 
+ T nn e iks on7„. 



*.(».*)= .^'""l. r".'"""^"' (7) 



Here T nn are the reflection coefficients, and T nm are the transmission coefficients. 
Substituting ip n into the matching conditions (jSJ) or (jSJ) we derive the scattering matrix 
that can be expressed via the coupling function 

/ 6\{a)-\ e x {a)B 2 {a) 6 1 {a)d 3 {a) \ 

S(a) = {-l) j - l 2 \ 9 1 (a)9 2 (a) 6 2 2 {a) - \ 6 2 {a)6 3 {a) 

\e 1 (a)e 3 (a) 6 2 (a)9 3 (a) 6 2 3 {a) - \ ) 

if a G Sj for j = 1,2. Note that the scattering matrix does not depend on k. 

Next let us look in detail at stationary scattering for the Hamiltonians H e (a, Q) 
and Ho and find the limit as e — > of the scattering amplitude. Consider the incoming 
monochromatic wave e~ lks coming from infinity along the edge 73. 

We shall seek a positive-energy solution of the problem 

-y" + aQ £ y = k 2 y on T \ {a}, VuJi (a) = y^a) = y^a), ^2—( a ) = Q ( 8 ) 

that coincides on Y £ \ Q E with ip 3 given by ((71). Let u a = u a (£, x), v a = v a (£, x) and 
w a = w a (£, k) be a linearly independent system of solutions of the problem 

- g " + aQg = K 2 g on Q \ {b}, 9ui {b) = g„M = gua (b), ^2^^ = ° ^ 

3=1 3 

on the graph Q. Clearly, the functions u a ((x — a)/e,ek), v a = v a ((x — a)/e,ek) and 
w a = w a ((x — a)/e,ek) form a linearly independent system of solutions of the problem 
(JBJ) on Q £ . Hence the desired solution on Q £ can be represented as a linear combination 
of this functions with the coefficients A, B and C respectively. Write x = ek. Since 
the potential aQ £ has discontinuities at the points a\, a\ and a|, we demand that the 
conditions 

y^) = y^), 0(o|) + ^(o|) =0, 3 = 1, 2, 3 (10) 

hold. On substituting our solution into f fTOj) we obtain the linear system Mc = r for 
the vector of unknown coefficients c(a,x) = (T 3 i,T 32 ,T 33 ,A,B,C) T . Here r(x) = 
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(0,0,0,0,e-" 



M(a, x) 



ixe 
( 



-ix\~V 



and 



V 



-ixe 










_ e ix 

—ixe 1 










_ e iK 

— ixe 1 



u a {ai,x) 

fe(«i>* 

u a (a 2 ,x) 

du n 



dui2 



a 2 , x 



u a {a 3} x) 

fe(«3,*) 



v a (ai,xj 

fe(«i>*) 

v a (a 2 ,x) 

fete.*) 

v Q (a 3 ,x) 

fete,*) 



w a (a 1 , 


x) 


\ 


dw a I 


x) 




w a (a 2 , 


x) 




dw a ( „ 


x) 




u> Q (a 3 , 


x) 




dlU Q ( 


x) 


/ 



Clearly, the functions u a (-,0), v a (-,0) and w a (-,0) form a linearly independent system 
of solutions of the problem ([2]) , ([3]) • In what follows we shall omit the second argument 
of these functions and keep in mind that it equals zero. We introduce the determinants 



hAa] 



ta -M fete) frte) 



fete) fete) £K) 

dn a 

, which is just the determinant h±- n with the 



^(a 3 ) 



d Sr(«s) 



dw, 

dui2 
dw, 



t(«3) 



u a (ai) u a (ai) iu a (ai) 

/in(«) = M a (a 2 ) v a (a 2 ) w a (a 2 ) 
u a (a 3 ) v a (a 3 ) w a (a 3 ) 

Let us also consider the function h 

m th -Yow replaced by the m th -row of h n for n = 0, 1; m = 1, 2, 3. 

Let us denote by A(x, a) the determinant of the matrix M(x, ct). It admits the 
asymptotic expansion 

A(x, a) = h(a) + ix{3/ii(a) - H (a)} + x 2 {-9/2 ^(a) + 3H (a) - H x (a)} + 0(x 3 ) 



as x — > 0, where i^, 
derive 



a 



E 



=1 /i„ m (a) for n 



0, 1. Employing Cramer's rule, we 



T 31 (x,a) 



A(x,a) 



u a {ai) 

fe(«0 

fete) 



Ua(fll) 

fete) 
fete) 



dw 



A(x,a) 



d'Wg 

dcU2 



T 32 (x,a) 



2i>c(l+ix) 
A{x,a) 



dUg 



Ol 



<foq 



w Q (a 2 ) 
fete) 



Ol 



w a (a 2 ) 
fete) 



loij 
(02) 

W a (Ol) 

fete) 
u a (a 2 ) 

fete) 

w a (a 2 ) 

fete) 



u a (oi) w a (ai) 



fete) 

v a {a 2 ) 



*£( 0l ) 



du>i 



w a (a 2 ) 



'M) 



0( 



X 



2;^ 
A(*,a) 



dtJ2 
Wa(Ol) 

w a (a 2 ) 

du 



u a (oi) iy a (oi) 



d^2 



O2J 



^33(x,a) 



-1 



(A(x,a)) 
-x 2 {#!(a 



/ii(a) + ix{H (a) 



v a {a 2 ) 
fe(*0 

2/i 03 (a) 



w a (a 2 

dw 



dui2 



a 2 



h{a)} 



+ flo(a) - 2(Ma) + M«)) - l/2/n(a)} 



0(x 3 



0(x 3 



as x —?• 0. Similar asymptotics can be derived without difficulty in the case, when 
the wave package comes from infinity along 71 or 7 2 . We shall study the asymptotic 
behavior as x — > of the scattering coefficients. Three cases are to be distinguished: 
A(x, a) = 0(1), A(x, a) = O(x) and A(x, a) = 0(x 2 ) as x -» 0. 
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Lemma 2.2 The set of roots of the equation h\(z) = coincides with the resonant set, 
i.e., hi is a characteristic determinant of the eigenvalue problem (0)-0). 

Proof. It is easily seen that the system u a , v a and w a of solutions of the problem (J9]) 
may be chosen so that 

«"<»> = "■ ;> = °- %® = l - 

Clearly, w Q vanishes on tui, and f a vanishes on u 2 . Let us consider the linear combination 

CLV a . . CLW a . . UMq, . . UW a . . fllto, . . Ufa . . . . 

# = 3 — (aiH — (a 2 )w Q + -, — (a 2 )- — {ai)v a - - — (a 2 )- — {ai)w a . (12) 
dui du 2 du 2 dui du 2 duj\ 

By construction ^f-(ai) = ^f"( a 2) = and ^-(03) = /11(a). Hence if h\(a) = and g 
is nontrivial on f2, then g is an eigenfunction of the problem (J5J), corresponding to a, 
and therefore a £ S. 

Note that at least one of the values f^(%), ^r( a i) or ifr( a j) * s different from 
zero for j = 1,2,3. Now consider the exceptional case, when g is trivial, i.e., all 
coefficients in (112j) equal zero. We show then that both values ^^(02) and ^r( a i) 
are zero. Conversely, suppose that f^(ai) 7^ 0. Since ^(01)^7(02) = 0, one 
obtains f^;(a 2 ) = fj;(a 2 ) = 7^(02) = 0, a contradiction. Similarly, f^(a 2 ) = 0. 
It follows that ^(aj = ^(a 2 ) = and f^(ai) = f£(a 2 ) = 0. Thus the function 
9i = fe(«sK " fe(os)«« Satisfies gL( ai ) = fe(a 2 ) = |i(a 3 ) = 0. If 9l is nontrivial 
on Q, then it is an eigenfunction of the problem (J2J), corresponding to a. If not, then 
u a would be a desired eigenfunction. 

What is left is to show that a £ £ implies h\(a) = 0. If we choose the new base 
u a , v a and w a such that u a is an eigenfunction of the problem (J2])-(|l]) corresponding to 
a, then the first column of the determinant is zero and the assertion follows. □ 

Lemma 2.3 The set of roots of the equation Hq(z) = that belongs to the resonant set 
coincides with E 2 . 

Proof. Suppose, contrary to our claim, that there exists a £ Si such that 
Ho(a) = 0. We choose the new linearly independent system u a , v a and w a such 
that u a is an eigenfunction of the problem ©-(jl]) corresponding to a. Write V = 
(feK^M'feM) and W = (fe(a0, £j(a 2 ), ^(a 3 )). These vectors are 
linearly independent. Indeed, in the opposite case there exists a nonzero number A 
such that V = XW. Therefore v a — \w a is an eigenfunction of the problem (J2J) (|4j) 
corresponding to a that is linearly independent with u a , a contradiction. 
In view of the Lagrange identities 

u a {a 1 )—^{a 1 ) + u a (a 2 )—^(a 2 ) + u a (a 3 )—^(a 3 ) = 0, 
du>x du 2 duj3 

u a {ai)—-^{a 1 ) + u a (a 2 )—-^(a 2 ) + u a (a 3 )— -(a 3 ) = 

CLUi CLU 2 UW3 
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the vector U - 
identity H (a) 



[U r 



(ai) , u a (a 2 ) , u a (a 3 )) is orthogonal to both vectors V and W. The 



can be written as 



u a {a lj 



du>2 
dw a 

dw2 



,02 



M &(«3) 
fe(*0 



+ w a (a 2 , 



(fog 



fe(«0 



fe(«s) fe(ai 



+^a(a 3 ) 



fe(«l) 



fe(«2) 



<*1, 



dw a 
du>2 



2 



0. 



It follows that [/ is orthogonal to the vector product V x W, hence that U is a zero 
vector, i.e., u a (oi) = ^ a (a 2 ) = ^(03) = 0, and finally that u a is a zero function, which 
is impossible. 

Let a belong to E 2 . We select a linearly independent system u a , v a and w a such 
that « a and v a form a base in the corresponding eigenspace. Obviously, the determinant 
Ho(a) vanishes, and the proof is complete. □ 



Remark 2.4 The roots of the equation hx(z) = as well as Hq(z) = do not depend 
on the choice of the linearly independent system of solutions of the problem (TJJ). 

Lemma 2.5 The function Hi is different from zero on E 2 . 

Proof. On the contrary, suppose that there exists a e £ 2 such that Hi(a) = 0. Let 
us select the system u a , v a and w a of linearly independent solutions of the problem (jHD 
such that u a and v a are eigenf unctions of the problem (J2J) (SD corresponding to a. Set 
U = (u a (ai), w a (o 2 ), u a (a 3 )) and V = (v a (ai), v a (a 2 ), v a (a 3 )). Note that the vectors U 
and V are linearly independent, since in the opposite case the functions u a and v a would 
be linearly dependent. 

From the Lagrange identities 



/ \ OjW a , . . . ctw a , . , . aw a , . 

Ua{ai)- — (01J + u a (a 2 )- — (a 2 ) + Ua(a 3 J- — (03) 
awi aa; 2 aa;3 

du> a , s , „ dw„ , „ , , dw; r 







v. 



(ai) j — («i) + Va(fl2) -3 — (02) + v a (a 3 )- — (a 3 ) = 



da;i 



da; 



it follows that the vector W 



' dw a 



Ol, 



dui 



«2 



The identity Hi(a) = can be written in the form 



duj 3 

a 3 )) is orthogonal to U and V. 



dw a 
duii 



oi, 



u Q (a 2 ) u a (a 3 ) 
v a (a 2 ) w a (a 3 ) 



diy 



(a 2 ) 



Ma(03) M Q (ai) 



v a (a 3 ) v a (ai 
dw n 



+ 



dw 3 



(as) 



u a (ai) u a {a 2 ) 
v a (a x ) v a (a 2 ) 



0. 



Consequently, W is orthogonal to the vector product U x V. Analysis similar to that 
in the proof of the previous lemma gives the assertion of the lemma. □ 

We are now in a position to prove the main result. Without loss of generality we 
establish the convergence of the transmission coefficient T 31 (ek,a) as e — > 0. Other 
coefficients may be handled in much the same way. 
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Proof of Theorem 12.11 The non-resonant case. Since a is not a resonant coupling 
constant, h\{a) ^ 0. From (fill) it immediately follows that T 31 (e/c, a) = 0(ek) as e — > 0. 
The case of simple resonance. If a £ Si, then hi(a) = and Ho(a) ^ in light of 
Lemmas [Q [231 By fTTT]) 



T 3 i(£fc,a) 



#o(a) 



w Q (oi) u a (ai) w a (ai) 



&« 



fe(«o 



&« 



feM fe(« 2 ) fe(a 2 ) 



O(efc), 



e -+0. 



We continue by choosing the new linearly independent system u Q , i> a and «;„ such that 
u a is an eigenfunction of the problem ([2])-(jl|) corresponding to a satisfying the condition 
u 2 a (ai) + u 2 a (a 2 ) + u 2 a (a 3 ) = 1 and 4^(a 2 ) = 0. Observe that 



T 31 (0,«) = 2fl 1 (a)^(a 1 )^(a 2 )tf - 1 i 



^^ d*™ " u ^ (13) 

Now let us express T3i(0, a) via the coupling function. Taking into account the Lagrange 
identities we see that 



u a {ai) 
u a (a 3 ) 



fefaO 



d^3 



II, 



a 2 )__fe(«3)Sf(«i)-^(«i)&(«3) 



fe(ai)' 



Therefore 
Ooi(a) = -w a (ai' 



dt> Q 



a 3 



dw c 



M 



« 



M a (a 3 ) 



dVct( n 



dw 2 



(02) 



\ 3 ( a 3j^ l°l, 



3 (a) 
(ai)- — (a 3 ) 



/102(a) 

L aa; 3 acui dw\ du 3 

By construction /i 03 (a) = u 2 a (a 3 )^(a 1 )^(a 2 )/ 
^02 (a) + /103(a)- Combining f[T4"|) and (fT5l) yields 



<M£(oi)SM 



03(a) 



(14) 



(15) 



Recall that Hq(cy) = /101(a) + 



#n 



feMfete) 



{w^(ai)+w^(a 2 )+^(a 3 )} 



da; 



*(*!)&(*.) 



3 (a) ^ - - ' - - - "'-- 03 ( Q ) 

Then by (fl3"l) 

T3 1 (O,a) = 20 1 (a)03(a) = r 31 (a). 

T/ie case 0/ double resonance. Finally suppose that a £ E 2 . By Lemmas 12.21 12.31 
and 12.51 we have h\{a) = Hq(q) = and Hi(pt) 7^ 0, thus 

Wa(oi) u a (ai) iw a (oi) 

fe(«i) &(«0 Sf(°0 + °( £fc )> £ ^° 

^0(02) Uq,(o 2 ) w a (a 2 ) 

by (fTTT) . Let us choose the new base u a , v a and w a such that w a and v a are eigenf unctions 
of the problem ([2]) corresponding to a satisfying 6 2 (a)+6 2 (a)+6 2 (a) = 1. Consequently, 

dw r 



T 3 i(eA;,a) 



#i(a) 



T 31 (0,a) = -2^ 3 (a)-^(a 1 )i7r 1 

OCJi 



a 
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Since the Lagrange identities hold it follows that 

fe("l) fefo) ^M 

0i (a) 5 2 (a) 5 3 (a) ' 

We thus get 

Ain(a) = e'i(a)- — (oi), /112(a) = 6> 2 (a)- — (a 2 ) = T-p-r-: — (oi), 

dui du2 "lKp.) u^i 

dw; a el{a)dw a 

Ma)=es(a)^(a 8 ) = ^ ? ^(a 1 ). 

By recalling Hx{a) = /in (a) + /112(a) + /113(a), one obtains 

HM) = ^vK(«) + ^(«) + ^(«)} = ^ff , 
hence 

T 31 (0,a) = -29 l {a)9 z {a) = T 31 (a), 
and the proof is complete. □ 



3. An example 

We suppose that the shape of a short-range potential in an actual model can be 
approximately described as follows 

{7 on wi, s G [0,0.5), 
-7 on wi, s e (0.5,1], 
otherwise. 

We study the scattering properties of the potentials ae~ 2 Q((x — a)/e) on Y £ as e — > for 
these models. As shown before, the scattering amplitude for the Hamiltonians H e (a, Q) 
and Ho converges as e — > to that for the limiting Hamiltonian H(a, Q). Certainly, the 
resonant set and the coupling function to be found are specific to the given shape Q. 

Table 1. Resonant intensities and coupling function 
a 61(a) 2 (a) =03 (a) 



8.8104 -0.9992 0.0279 

28.5513 0.9999 0.0012 

59.5701 -0.9999 0.997 x 10~ 4 



Since the function Q vanishes outside the edge uj\, the resonant set £ coincides with 
the set Si of the simple eigenvalues of (E])-(j4j). In fact, all eigenfunctions are constant 
outside the edge U\ and the resonant set is the spectrum of the problem 

. dg dg 

-g + aQg = on uj 1 , - — (oi) = 3 — (6) = 0. 

dui dcji 
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Table 2. 


Scattering matrix for a = 8.8104 


0.9968 

-0.0558 
-0.0558 


-0.0558 -0.0558 
-0.9984 0.0016 
0.0016 -0.9984 


Table 3. 


Scattering matrix for a = 28.5513 


0.9996 
0.0024 
0.0024 


0.0024 0.0024 
-0.9999 0.288 x 10~ 5 
0.288 x 10- 5 -0.9999 


Table 4. 


Scattering matrix for a = 59.5701 


0.9996 
-0.0002 
-0.0002 


-0.0002 -0.0002 
-0.9999 0.199 x 10~ 7 
0.199 x 10~ 7 -0.9999 



Table [T] lists the first three positive resonant values of a (numerically computed 
using Maple) and the corresponding values of the coupling function 9(a). Tables [2J [3] 
and S] list the scattering matrixes for different resonant intensities. We note that the 
transmission coefficients decay very fast. 

The similar effects were observed, when the asymptotic behavior of the scattering 
coefficients for K £ (a,Q) was discovered directly to the given shape Q. Figure |2] shows 
the logarithm of the transmission probability \T 3 i(ek, a)\ 2 versus a for ek = 0.0001 and 



8 



at) 
o 




Figure 2. Logarithm of the transmission probability \T^,i(ek , a)\ 2 against a at 
ek = 0.0001 (solid curve) and 0.01 (dash curve) 
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ek = 0.01. From this figure we see that \T 3 i(ek, a)\ 2 vanishes as e — > except when a 
takes resonant values which correspond to the spikes of |T 31 (e/c, a)\ 2 . Observe also that 
the transmission coefficient T^i(ek,a) decays as a — > oo. 
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